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Abstract 

Background: Theileria and Babesia protozoan parasites are transmitted mainly by ticl< vectors. Tliese parasites cause 
lieavy economic losses to the live-stock industry, as well as affecting the health of wild animals in parasite-endemic 
areas. Identification of infectious agents in wild animals is not only crucial for species preservation, but also provides 
valuable information on parasite epidemiology. Here, we conducted a molecular surveillance study in Northwestern 
China to assess the prevalence of blood pathogens in cervids. 

Methods: PGR analysis and microscopic evaluation of blood smears to detect Theileria- and Babesia-related diseases 
in Cervidae were conducted, in which 22 blood samples from red deer (n = 22) in Qilian Mountain and 20 from sika 
deer (n = 20) in Long Mountain were collected and tested for the presence of Theileria and Babesia. The 18S rRNA 
gene was amplified, and selected polymerase chain reaction (PCR)-positive samples were sequenced for species 
identification. 

Results: PCR revealed that 9.1% of the Qilian Mountain samples and 20% of the Long Mountain samples were 
positive for Theileria uilenbergi; 90.09% of the Qilian Mountain samples (n = 22) were positive for T capreoli, but all of 
the Long Mountain samples (n = 20) were negative for T capreoli; no other Theileria or Babesia species were found. 
PCR showed that T uilenbergi and T capreoli were present in red deer in Qilian Mountain, while only T uilenbergi 
was found in sika Deer in Long Mountain. The 18S rRNA gene sequences were aligned against the corresponding 
GenBank sequences of known isolates of Theileria and Babesia and subjected to phylogenetic analysis. The 
phylogenetic tree showed that the newly isolated Theileria spp. could be classified as belonging to two clades: one 
group belonged to the same clade as T uilenbergi, the other to a clade containing T capreoli. 

Conclusions: Our results provide important data to increase understanding of the epidemiology of Cervidae 
theileriosis, and will assist with the implementation of measures to control theileriosis transmission to Cervidae 
and small ruminants in central China. 
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Background 

Theileria and Babesia, parasites that are mainly transmitted 
by tick vectors, cause heavy economic losses to live-stock 
and affect the health of wild animals where ever such para- 
sites are endemic. Theileria and the closely related genus, 
Babesia, exhibit complex life cycles that involve mamma- 
lian intermediate hosts and hard ticks as the definite hosts 
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[1]. Theileria is an obligatory intracellular parasite that in- 
fects leukocytes and erythrocytes in the intermediate host, 
unlike Babesia, which only infects erythrocytes. Among the 
Theileria and Babesia parasites that infect cervids, several 
Theileria spp. (e.g. T cervi, T capreoli, Theileria sp. OT3, 
Theileria sp. ZS OT4, and T ovis), and Babesia spp. (e.g. B, 
bigemina, B. bovis, B. capreoli, B. divergens, and B, odocoi- 
lei) are described as being moderately pathogenic or benign. 
Theileria and Babesia are cosmopolitan parasites [2] that 
have been detected in wild ruminants in many countries in- 
cluding Japan [3-5], South Korea [6], Brazil [7,8], the United 
States [9-12], Italy [13], eastern Austria [14], northern and 
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central Spain [15-17], and the United Kingdom [18]. In 
China, He et al, [19] first reported the existence of Theileria 
infection in sika deer in the Hubei province of central 
China. 

Although there are no published epidemiological data 
on the Cervidae of northwestern China, no suspected 
cases of theileriosis in small ruminants have been ob- 
served here. In the current study, we performed an epi- 
demiologic survey of Theileria and Babesia infections of 
Cervidae in northwestern China. Polymerase chain reac- 
tion (PCR) analysis and microscopic examination of 
blood smears were used for species identification and to 
elucidate the evolutionary relationships among the newly 
identified hemoparasite species based on the 18S rRNA 
gene sequences of known Theileria and Babesia species. 

Methods 

Sample collection 

The two regions investigated in northwestern China are 
located at latitudes 34° 44' to 73° 28' north and longitudes 
97° 20' to 106° 35' east. The study period was August 



2013. Twenty-two blood samples (n = 22) were collected 
randomly from wild red deer in Qilian Mountain, and 20 
blood samples (n = 20) were collected randomly from 
domesticated sika deer at a ranch in Long Mountain, 
northwestern China. Blood smears were prepared for both 
groups (Table 1). During the blood collection process, cases 
of suspected theileriosis and babesiosis were investigated. 
Theileriosis and/or babesiosis should be suspected in tick- 
infested animals with a fever, enlarged lymph nodes (only 
for theileriosis), anemia and jaundice, or hemoglobinuria 
(only for babesiosis). 

Microscopy of blood smears 

Blood smears were air-dried, fixed in methanol, stained 
with a 10% solution of Giemsa in phosphate-buffered sa- 
line (pH 7.2), and then subjected to microscopic analysis 
and photography. 

DNA extraction 

Genomic DNA from 42 whole blood samples was ex- 
tracted using a genomic DNA extraction kit (Qiagen, 



Table 1 Sequences of the oligonucleotide primers used in this study 

Pathogen Target gene Oligonucleotide sequences (5-3') 



Forward 



Reverse 



Final amplicon Reference 
size (bp) 



Hemoparasite 


18S rRNA 


Theileria spp. 


18S rRNA 


T. onnuloto 


30 kDa 


T. sergenti 


MPSP 


T. sinensis 


MPSP 


T. oriental is 


MPSP 


T. cervi (Type 


18S rRNA 


F) 




T. cervi (Type 


IBS rRNA 


G) 




T. luwenshuni 


18S rRNA 


T. uilenbergi 


18S rRNA 


T. ovis 


IBS rRNA 


T. lestoquardi 


18S rRNA 


T. copreoli 


18S rRNA 


Babesia spp. 


IBS rRNA 


B. odocoilei 


18S rRNA 


B. bovis 


18S rRNA 


B. ovata 


AMA-1 


B. bigemina 


IBS rRNA 


B. ovis 


18S rRNA 


B. motasi 


18S rRNA 



AACCTGGTOATCCTGCCAGT 
AGmCTGACCTATCAG 

GTAACCmAAAAACGT 

CACGCTATGTOTCCAAGAG 
CACTGCTATGTOTCCAAGAGATA^ 
CmGCCTAGGATACTOCT 
TOCCmGAGGGGT 

OTCCCGTOTGGAGG 

GGTAGGGTATOGCCTACTGA 

GGTAGGGTATOGCCTACCGG 

TCGA-GACCTOGGGT 

GTGCCGCAAGTGAGTCA 

gtcttgtaattggaatgatgg 

TGTOTGAATAOT(C/G)AGCATGGAA 



GATCCTOTGCAGGTOACCTAC 
TOCCmAAOTCCTO 

GTOCGAACATGGGm 

TGTGAGACTCAATGCGCCTA 
AATGCGCCTAAAGATAGTAGAAAAC 
ACGGCAAGTGGTGAGAACT 
AAGCCTATOCCGTACC 

CTOACCTATOCCGTGAG 

TCATCCGGATAATACAAGT 
ACACTCGGA AAATGCAAGCA 
TCCGGACATOTAAAACAAA 
GGACTGATGAGAAGACGATGAG 
ccaaagactttgatttctctc 

CGAOTCTCCmAAGTGATAAC 



CACCGTA^GAC^mGTCGACTGTCGG CCCGTAACGGACGAACCTOTCACGGG 
CTGTCGTACCGTOGTOAC CGCACGGACGGAGACCGA 
GATACGAGGCTGTCGGTAGC AGTATAGGTGAGCATCAGTG 



TGGCGGCGmATOGTCG 

TGGGCAGGACCTOGTOTOT 

AAGAAmCACCTCTGACAG 



CCACGOTGAAGCACAGGA 

CCGCGTAGCGCCGGCTAAATA 

GOTGC^mGTOCmG 



1750 
1098 

721 

875 
887 
776 
654 

654 

340 
340 
520 
785 
580 

950 
900 
541 
504 

1124 

549 

179 



Medlin et al., [20] 

Allosop et al., 
[21] 

d'Oliveria et al., 
[22] 

Liu et al., [23] 
Liu et al., [23] 
Ota et al., [24] 
Boes et al., [25] 

Boes et al., [25] 

Yin et al., [26] 

Yin et al., [26] 

Altay et al., [27] 

Aktas et al., [28] 

Slemenda et al., 
[29] 

Ramos et al., [10] 

Ramos et al., [10] 

Guido et al., [30] 

Sivakumar et al., 
[31] 

Guido et al., [30] 

Aktas et al., [32] 

Shayan et al., 
[33] 
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Hilden, Germany) according to the manufacturers in- 
structions. The DNA yields were determined using a 
NanoDrop ND-2000 Spectrophotometer (Nanodrop 
Technologies, Wilmington, DE, USA). 

Molecular detection of Theileria and Babesia using 
species-specific primers 

PGR was used to detect and differentiate Theileria and 
Babesia using species-specific primers [10,20-32], the 
details of which are shown in Table 2. PGR reactions 
were performed in an automatic DNA thermocycler 
(Bio-Rad, Hercules, GA, USA) and PGR products were 
separated by 1.5% agarose gel electrophoresis to assess 
the presence of specific bands indicative of Theileria 
spp. and Babesia spp. 

PGR amplification of the 18S rRNA gene of Theileria 

PGR was used to amplify protozoan 18S rRNA gene se- 
quences with primers A/B as described by Medlin et al 
[20]. The DNA fragments generated were ligated into 
pGEM T easy vectors (Invitrogen, Garlsbad, GA, USA) and 
transformed into competent JM109 cells (Takara Bio Inc., 
Shiga, Japan). The transformed bacteria were plated onto 
selective LB medium, and at least three positive clones 
identified by the 989/990 primers specific to Theileria spp. 
[21] were sequenced by the GenScript Gorporation (Piscat- 
away, NJ, USA). Representative sequences of the newly 
identified Theileria and Babesia IBS rRNA genes from this 
study were deposited in the GenBank database of the 
National Genter for Biotechnology Information (NGBI) 
(http://www.ncbi.nlm.nih.gov/genbank/). 

Sequence alignments and phylogenetic analysis 

Gompilation, editing, and assembly of the multiple se- 
quences generated from each template were performed 
using the EditSeq and SeqMan algorithms of the Laser- 
gene software package for Windows (DNASTAR, Madi- 
son, WI, USA). Sequence alignment and phylogenetic 
analysis performed by the MegAlign component of the 
Lasergene program ver. 4.01 (DNASTAR) were used to 
perform multiple sequence alignments with the Glus- 
talW algorithm and phylogenetic analysis by the 
neighbor- joining method. A phylogenetic tree was con- 
structed (Figure 1) based on the Theileria and Babesia 
18S RNA gene sequences determined in the present 
study, other sequences from our laboratory, and those 



obtained from GenBank under the following accession 
numbers: KJ188207-KJ188232, AY726011, DQ866842, 
AY421708, AY661512, AY661515, EU274472, EU277003, 
AY260171, AY260172, FJ603460, EU083800, FJ426369, 
M64243, AY262116, AY262118, AY262119, JX469515, 
AY262120, JF719835, FJ595120, FJ599640, AB012201, 
AB012194, HQ188406, AF529272, Z15105, AY260176, 
AY260177, AY260178, AY260179, AY260180, AB576641, 
AB012198, GQ304524, AY260180, AY081192, HQ264111, 
and HQ264112. In the phylogenetic tree, the length of each 
branch pair represents the distance between sequence pairs, 
while the units on the horizontal axis indicate the number 
of substitution events (Figure 1). 

Ethical approval and consent 

This molecular approach was then employed to identify 
Theileria parasites of northwestern Ghinese cervidae at 
Lanzhou Veterinary Research Institute, Ghinese Academy 
of Agricultural Science. Use of these clinical samples was 
approved by the National Review Board (Ghina). 

Results 

Microscopic examination of thin blood smears 

Piroplasm infections in the blood smears were observed 
microscopically in samples from Qilian Mountain (72.7%, 
16/22) and Long Mountain (10%, 2/20); all infections 
exhibited low parasitemia levels (0.01-2.5%) (Table 2). 
The piroplasms were polymorphous, and most presented 
as rod-like, needle-like, pear-shaped, or spherical forms 
(Figure 1). Gases of theileriosis and babesiosis were not seen 
in the Qilian Mountain and Long Mountain regions of 
northwestern Ghina. 

PGR detection of Theileria and Babesia using species- 
specific primer sets 

All 42 blood samples were negative for Babesia spp. 
using Babesia genus and Babesia species-specific primer 
sets; 24 of them (n = 42) were positive for Theileria spp. 
using Theileria genus and Theileria species-specific pri- 
mer sets. PGR using twenty sets of primers revealed that 
only two Theileria species, T, capreoli and T, uilenbergi, 
were detectable in the blood samples. The PGR-positive 
rates for T. uilenbergi in the blood samples from red 
deer in Qilian Mountain and sika deer in Long Moun- 
tain were 9.1% (2/22), and 20% (4/20), respectively. 
While the PGR-positive rates for T. capreoli in blood 



Table 2 Prevalence of Theileria and Babesia in red deer and sika deer in northwestern China, 2013 



Collection 


Host 




Prevalence of Theileria spp. and Babesia spp. in blood samples 


sites 




ME* 


T. capreoli T. uilenbergi Babesia spp. 


Qilian Mountain 


Red deer 


72.7% (16/22) 


90.9% (20/22) 9.1 % (2/22) 0% (0/22) 


Long Mountain 


Sil<a deer 


10% (2/20) 


0% (0/20) 20% (4/20) 0% (0/20) 



*ME, microscopic examination of blood smears. 
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Figure 1 Piroplasms in blood smears from Cervidae. 



samples from red deer in Qilian Mountain were 90.9% 
(20/22), T. capreoli was not detected in the blood sam- 
ples from sika deer in Long Mountain (Table 2). These 
results confirm that all of the samples that were blood 
smear-positive were also PCR-positive for Theileria par- 
asites. Of the 24 samples for which no piroplasms were 
observed by microscopy, we found that one was 
r. capreoli- and T, uilenbergi-positive, three were T. 
capreoli'positive, and two were T. uilenbergi-positive by 
PGR. Compared with microscopy of blood smears, the 
PGR method had higher sensitivity and was recom- 
mended to be utilized in the diagnosis of piroplasmosis. 

18S rRNA amplification 

The 18S rRNA gene sequences were determined for the 
Theileria spp. isolates from Qilian Mountain and Long 
Mountain. The length of the 18S rDNA gene was 
1,751 bp (r. capreoli) and 1,745 bp {T, uilenbergi) using 
the A/B primers, which are specific for hemoparasites. 
The 18S rRNA gene sequences obtained from red deer 
and sika deer in northwestern China have been depos- 
ited in GenBank under the following accession numbers: 
KJ188207-KJ188220 for the T. capreoli isolates from 
Qilian Mountain red deer; KJ188227-KJ188232 for the 
T, uilenbergi isolates from Qilian Mountain red deer; 
and, KJ188221-KJ188226 for T. uilenbergi isolates from 
Long Mountain sika deer. 

Sequence alignments and phylogenetic analysis 

The Basic Local Alignment Search Tool (BLAST) at NCBI 
was used to analyze the sequences obtained in our study. 
The results showed that there were two Theileria se- 
quences from Cervidae in Qilian Mountain; one shared 



99-100% homology with Theileria sp. 385/02, T capreoli 
isolate BAB 11 58, and Theileria sp. CNY3B, while another 
shared 100% homology with the following isolates: T uilen- 
bergi Longde, T. uilenbergi Yongjing, and T. uilenbergi 
Zhangchuan. However, only one species of Theileria was 
found in Long Mountain Cervidae, and its 18S rRNA gene 
sequence shared 100% homology with the T. uilenbergi 
Longde isolate, the T. uilenbergi Yongjing isolate, and the 
T. uilenbergi Zhangchuan isolate. Sequence alignments 
showed that the newly identified Theileria spp. sequences 
with accession numbers KJ188207-KJ188220 were highly 
homologous to each other with identity values of 100%, 
while the other set of newly identified Theileria sequences 
(KJ188221-KJ188232) also were highly homologous to each 
other with identity values of 100%. Consistent with the 
BLAST results, the phylogenetic tree constructed in this 
study showed that one newly sequenced Theileria isolate 
appeared in the same clade as T. capreoli, while another ap- 
peared in the same clade as T. uilenbergi (Figure 2); both of 
the new 18S rDNA gene sequences show obvious differ- 
ences to those from Babesia (such as those from B. ovis, 
B. bovis, B. motasi, B. ovata, and B, divergens) or other Thei- 
leria sequences (such as those from T ovis, T luwenshunli, 
T sinensis y T buffei/T sergenti, T equi, and T annulata). 

Discussion 

To the best of our knowledge, this study is the first to 
report on cervine theileriosis in northwestern China. 
Molecular screening for Theileria and Babesia in red 
deer and sika deer in northwestern China showed that 
the most prevalent species of Theileria in the deer were 
T, capreoli and T, uilernbergi; however, no Babesia para- 
sites were found. The prevalence of T, capreoli was very 
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J AY260172XOVI5 Turkey 
f FJ603460 Tovis China 
AY260171 T.ovis Sudan 

EAY2621 18 T hiwenshum Ningxia 
JX469515 T. hiyfenshwn Songxian 
AY2621 19 T. bxyfenshum Qinghai 

Y}595m Theilena sp Hani 

Yimm Theilena sp Han2 

ABO 12201 Theileria sp. CC3A 

P-r~ AY661515X sergentx Japan- Iked 
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I 1 EU274472 T. sinensis Lint an 

EU277003 T. swensis Weiyuan 
[— AB012194 Theilena sp CNYIA Japan 
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M64243 T. annulala 
EU083800X annulata Ningxia 
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KJ1S8221 Theileria sp. Long Mountain 1 
KJ188222 Theileria xp LougMouutaui 2 
KJ188223 Theileria sp. LongMountain 3 
IF7 19835 T. uilenbergi Yongjing 
KJi8822" Theileria sp. QiUaii Mountain 1 
I KJi88228 Theileria sp. Qilian Mountain 2 

[ KJ188229 Theileria sp. Qihaii Mountain 3 

AY262 1 20 T. uilenbergi Longde 

Z15105 T.equi South Afnca 

I AY2601795. woto57 Netherlands 
AY260180 B. motasi Netherlands Tex e 
kYm\92 B.ovata Korea 
AY260\76 B.crassa Iran 
AY260177 5. d crassa GU184Turkey 
AB 576641 5. wicwfit Japan 
GQ304524 5. divervens France 
AY2eom B.ovis Turkey 



-I — HQ2641115. bovts USA 
EQ2eAU2B.bovis USA 



16 14 12 10 8 6 4 2 

Nucleotide Substitutions (xlOO) 

Figure 2 Phylogenetic tree of Theileria and Babesia spp. based on IBS rRNA gene sequences. 



high in Qilian Mountain deer, unlike T. uilenbergU the 
prevalence of which was relatively low in Qilian Moun- 
tain. In contrast, the prevalence of T, uilenbergi was 
relatively low in Long Mountain, and no T, capreoli in- 
fections were found at this location. In both of the areas 
investigated, no Babesia spp. (including B, odocoilei, 
B. bovis, B. ovata, B. bigemina, B. ovis, and B, motasi) or 
infections with other Theileria spp. (including T, cervix 
T, luwenshuniy T, ovis, T, sergenti/T, buffeli/T, orientalise 
T. annulata, and T. sinensis) were detected using genus- 
and species-specific primers. It has been speculated 
that in areas with high infection rates, there is an in- 
verse relationship between age and resistance to infection, 
where fawns gradually acquire immunity without showing 
clinical symptoms, and immunity is maintained by repeated 
challenges with the parasites. Consequently, a persistent 
parasite reservoir is established in the wild ruminants [16]. 
However, in these hosts, stressors like high parasitemia. 



poor nutrition, high population density, harsh weather 
conditions, or handling (e.g. translocation) can lead to 
symptomatic piroplasmosis, a cause of severe disease 
and death [15]. 

The diagnosis of parasitic diseases and detection of para- 
site species is often difficult at the carrier stage and in 
animals with mixed infections. Nevertheless, molecular 
diagnostic methods can enable the direct, specific, and sen- 
sitive detection of parasite species [27,34]. PCR-based mo- 
lecular techniques allow the detection and discrimination 
of these parasites at low parasitemias and in animals with 
mixed infections [33,34]. The PGR methods used in this 
survey can detect cervine Theileria and Babesia directly in 
the same samples based on the 18S rRNA gene sequence 
using Theileria and Babesia specific primer sets. In the 
present study, the prevalence of Theileria infections in red 
deer and silk deer, which were detected by PGR, was 
statistically higher than that observed via microscopic 
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examination. These results were expected because carrier 
animals often have low parasitemias, which are dift'icult to 
detect microscopically. 

The maximum parsimony tree of the 18S rRNA 
genes indicated that the cervid infections detected in 
this study could be divided into two groups within the Thei- 
leria species clade. Group 1 sequences fell within the 
low-pathogenicity clade containing T, capreoli BAB1158 
(AY7266011: isolated from European roe deer in Spain), 
Theileria spp. 3185/02 (AY421708: isolated from red deer 
in Spain), and Theileria spp. 3185/02 (DQ866842: isolated 
from roe deer from the Basque Country of Spain), and were 
also present in roe deer and red deer from northern Spain 
[16], and in fallow deer from Italy [13], but absent from 
chamois. In contrast, group 2 Theileria spp. fell into the 
same clade as the high-pathogenicity T uilenhergi Longde 
isolate (AY262120), as did the T. uilenhergi Lintan isolate 
(AY262116) and the T. uilenhergi Yongjing isolate 
(JF719835), which are all widely present in sheep and goats 
and transmitted by Haemaphysalis qinghaiensis and H. 
longicornis in northern and northwestern China; their se- 
quences were almost identical to each other with high boot- 
strap values (Figure 1). Clustering of the 18S rRNA 
sequences into two groups indicates the heterogeneity of 
the 18S rRNA genes of cervine Theileria spp. in China. 
Therefore, the two Theileria spp. newly identified in Cervi- 
dae in northwestern China should be classified as T 
capreoli and T uilenhergi. The 18S rRNA gene sequence of 
T capreoli was most closely related to Theileria CYN3B of 
Japan and T ovis, which belong to an evolutionary distinct 
lineage of non-lymphoproliferative Theileria species. This 
species was clearly divergent from Bahesia and a second 
lineage of lymphoproliferative Theileria species that in- 
cluded T annulata and T parva. 

Ideally, tick transmission studies on ticks that are asso- 
ciated with infected herds should be conducted to iden- 
tibf the species involved in transmitting theileriosis to 
Cervidae. In this study, no ticks were collected from the 
deer or the surrounding environments of the two areas 
investigated because the sampling was carried out in Au- 
gust, a time of year when the ticks are inactive. Yin et al 
[35] and Li et al [36] identified H, qinghaiensis and H, 
longicornis as vectors for T, uilenhergi in China. H, qin- 
ghaiensis is a triphasic tick common in Northwestern 
China and H, longicornis is a ubiquitous triphasic tick 
found in most parts of China. Therefore, it can be spec- 
ulated that H, qinghaiensis and H, longicornis may play 
an important role as natural vectors of T, uilenhergi in 
northwestern China. Galuppi et al [13] detected T, 
capreoli DNA in Ixodes ricinus, but whether other vec- 
tors transmit T. capreoli remains unclear. At present, 
identification of the full range of tick species that can 
transmit T. capreoli to Cervidae awaits further 
investigation. 



Conclusions 

Our results provide important data that should increase 
understanding of the epidemiology of cervine theilerio- 
sis, and assist with the implementation of measures to 
control theileriosis transmission to Cervidae and small 
ruminants in northwestern China. Clarification of the 
role that Cervidae might have as reservoir hosts for 
maintaining T. capreoli and T. uilenhergi is critical to 
understanding whether Cervidae contribute to the 
spread of ruminant theileriosis in China. 
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